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We perform a search for the exclusive radiative decay B° — > <fry, which is dominated by bd 
annihilation, in a sample of 124 million BB events collected with the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage ring at SLAC. No significant signal is seen. We set an upper limit 
on the branching fraction of B(B° — > (fry) < 8.5 x 10 -7 at the 90% confidence level. 

PACS numbers: 13.20.Ho 



Within the Standard Model (SM) the rare decay B° -> 
dry proceeds through a penguin annihilation process. No 
process of this kind has yet been observed. The largest 
short-distance contribution to the SM amplitude is illus- 
trated in Fig. ^ The coupling of the top quark within 
the loop to the d quark leads to a dependence of the 
amplitude on the Cabbibo-Kobayashi-Maskawa matrix 
element V t d Q that suppresses the decay rate. Helic- 
ity suppression and the smallness of the decay constants 
/b and also lead to a small decay rate 2]. A recent 
QCD factorization calculation predicts a SM branching 
fraction of 3.6 x 10~ 12 2]. However, contributions to 
the B — ► </>7 amplitude are possible from physics be- 
yond the SM where new heavy particles enter the loop. 
For example, some models of supersymmetry that vio- 
late i?-parity predict an enhancement of the B° — > dry 
branching fraction by up to four orders of magnitude [2| . 
A prior experiment has bounded the branching fraction 
to be B(B° -> (jry) < 3.3 x 1(T 6 at the 90% confidence 
level (C.L.) ||. 

This analysis uses data collected with the BABAR de- 
tector at the PEP-II asymmetric-energy e + e~ storage 
ring at SLAC. The data sample consists of 124.1 ± 1.4 
million BB events, corresponding to an integrated lumi- 
nosity of 113 flV 1 on the T(4S) resonance, which has a 
mass of 10.58 GeV/c 2 . 

The BABAR detector is described in detail in Ref. 
Q- Charged particle trajectories are measured by the 
combination of a five-layer silicon vertex tracker (SVT) 
and a 40-layer drift chamber (DCH), which are embed- 
ded in the 1.5 T magnetic field of a solenoid. Pho- 
tons are detected in a CsI(Tl) crystal electromagnetic 
calorimeter (EMC) , with an energy resolution of <je /E = 
0.023( J B/GeV)- 1 / 4 © 0.019. Charged hadron identi- 
fication is performed by the combination of energy- 
loss (dE/dx) information from the SVT and DCH and 
measurements from a ring-imaging Cherenkov detector 
(DIRC). The segmented flux return (IFR) of the magnet 
is instrumented with resistive plate chambers to identify 
muons. 

Monte Carlo (MC) simulations of the response of the 
BABAR detector, based on GEANT4 0, are used to op- 
timize the selection criteria and determine the signal ef- 
ficiencies. These simulations take into account the vari- 
ations of detector conditions and beam induced back- 




FIG. 1: One of the leading order Feynman diagrams con- 
tributing to the decay B° — > (fry in the Standard Model. 

grounds during the data-taking period. 

The first stage of the analysis is to identify d> mesons 
and high-energy photons, which can be combined to form 
B° candidates [6J . The (j> is reconstructed in the decay to 
K + K~, which corresponds to (49.2^)% of the total cj) 
decay rate 0. We require that the charged tracks used 
in reconstructing the </> candidates have associated hits 
in both the SVT and DCH and have a transverse mo- 
mentum greater than 0.1 GcV/c. Tracks compatible with 
being kaons are identified via an algorithm that com- 
bines the candidate track's measured dE/dx, Cherenkov 
angle and number of Cherenkov photons. Individual 
kaons produced by signal <fi meson decays are identified 
with 80% efficiency by this algorithm, while the misiden- 
tification rate of pions as kaons is less than 2% over 
most of the relevant kaon momentum range. A d) can- 
didate is composed of two identified kaons of opposite 
charge that are consistent with originating from a com- 
mon vertex. We require d> candidates have a mass of 
1.011 < m K+K - < 1.029 GeV/c 2 (a full width 4.2 times 
the natural width); this criterion is optimized by a pro- 
cedure described below. 

A photon is identified as a shower in the EMC that 
is not associated with a reconstructed track. We re- 
move poorly reconstructed photons by rejecting show- 
ers in crystals with high noise rates or rejecting photons 
that are near the edge of the calorimeter's acceptance. 
Furthermore, showers that contain energy deposits mea- 
sured by crystals adjacent to an inefficient crystal are 
rejected. The shower profile is required to be compati- 
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FIG. 2: The <j> invariant mass distribution of on-resonance 
data. The <f> candidates are selected in the region indicated by 
the arrows. All selections have been applied except the jties, 
AE* , cj> mass, and neural net requirements. The distribution 
is dominated by combinatorial selections. 



FIG. 3: Neural network output for MC-simulated events with 
comparison to data control samples. The histograms have 
been normalized such that their total area is unity. The on- 
resonance sideband is defined as regions (I) and (II) in the 
text. The arrow indicates the lower limit imposed on the 
data. 



ble with a single photon to reject those that arise from 
7T° decays where the two photons enter the calorimeter in 
close proximity to one another; this condition also rejects 
showers generated by neutral hadrons, predominantly n 
and K®, that have a broader lateral profile than pho- 
ton showers. To suppress photons from ir° (ry) decays, 
the photon candidate is combined in turn with all other 
photons in the event with a laboratory energy greater 
than 50 (250) MeV. If any of the resulting invariant 77 
masses are within 20 (40) MeV/c 2 of the 7r° (77) mass, the 
candidate is vetoed. To further remove photons from tt° 
decays the shower is required to be isolated by at least 
25 cm from any other shower in the event. 

The photon and <f> meson candidates are combined to 
form B° meson candidates. We define AE* = E* B — 
^beami where E* B =E*^ + E* is the center-of-mass (CM) 
energy of the B° meson candidate and E^eam * s the 
CM energy of each beam. The signal AE* distribu- 
tion is peaked at zero with a resolution of approximately 
50 MeV; there is a negative tail in AE* because of the 
asymmetric E* resolution that is the result of energy 
leakage from the EMC. We also define the beam-energy- 
substituted mass m-Es = V^bcam — Pb ' w here p B is 
the CM three-momentum of the B° candidate. The sig- 
nal ?7Zes distribution peaks at the mass of the B meson, 
ms = 5.279 GeV/c 2 7], and has a resolution of 3 MeV/c 2 , 
which is dominated by the spread in E^ cam . The signal 
region in which we search for B — > <jyy events is de- 
fined as -0.2 < AE* < 0.1 GeV and 5.27 < m ES < 
5.29 GeV/c 2 . We define three sideband regions for esti- 
mating backgrounds: (I) 5.10 < mss < 5.29 GeV/c 2 and 
0.1 < AE* < 0.5 GeV, (II) 5.10 < to es < 5.29 GeV/c 2 
and -0.5 < AE* < -0.2 GeV, and (III) 5.10 < m ES < 
5.27 GeV/c 2 and -0.2 < AE* < 0.1 GeV. At this 



stage in the analysis all candidates with —0.5 < AE* < 
0.5 GeV and 5.10 < to E s < 5.29 GeV/c 2 are retained. 

The background comes predominantly from random 
combinations of real cf> mesons (Fig. [2Jl and high-energy 
photons produced in continuum u, d, s and c quark- 
antiquark events. The dominant sources of high-energy 
photons in continuum events are initial-state radiation 
(ISR) and ir°/r] — > 77 decays where the second photon is 
undetected, or the measured two-photon mass lies out- 
side the 7T° jr\ veto window. Since the continuum back- 
ground does not peak in AE* or m E s , its magnitude can 
be evaluated from parts of the mES and AE* data distri- 
butions in the sideband regions. There are also potential 
backgrounds from charmless B decays that peak in mES 
and AE* . The BB background is estimated from simu- 
lation. 

The ratio of the second-to-zeroth Fox-Wolfram mo- 
ments 0, for charged tracks in the event measured in 
the CM frame, is required to be less than 0.9 to reject 
some low multiplicity continuum final states. The com- 
binatorial background within the signal region is reduced 
further by combining 24 input variables that distinguish 
between qq continuum and BB events into one discrimi- 
nating variable via a neural network, as in Refs. [jl fiof. 
The neural network responds non-linearly to the input 
variables and exploits correlations among them To 
discriminate between jet-like continuum background and 
the more spherically-symmetric signal events, we include 
in the neural network the absolute value of the cosine of 
the angle between the high-energy photon and the CM 
thrust axis of the reconstructed particles in the r.o.e (rest 
of the event), | cos^|, and the energy distribution of all 
reconstructed particles in the r.o.e. binned into eighteen 
10° intervals around the photon direction. The distri- 
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bution of events in | cos 0*-, | is uniform for signal events 
but is strongly peaked toward 1 for continuum events. 
The energy flow is collinear with the photon direction in 
continuum events but is more isotropic in signal events. 
To discriminate against ISR continuum events we incor- 
porate in the neural network the ratio of the second-to- 
zeroth Fox- Wolfram moments computed in the photon 
recoil frame, in which ISR continuum events are more 
likely than signal events to have a back-to-back jet struc- 
ture. To provide further discrimination from continuum 
background, we include in the selection two angular vari- 
ables of the signal decay: | cos 0* B | , where 0* B is the an- 
gle between the B meson candidate's momentum and the 
beam axis in the CM frame, and | cos Oh | , where Oh is the 
angle between the flight directions of one of the daugh- 
ters of the tfi meson and the B meson candidate in the 
rest frame of the </> meson. The distribution of events of 
both of these variables is proportional to sin 2 0, where 
is B or Oh for signal events. The cos 6*^ distribu- 
tion has this form because a vector state (T(4S')) is de- 
caying to two pseudoscalars (BB), with the vector state 
having no helicity zero component. The \ cos 0h\ distri- 
bution has this form because a pseudoscalar state (B°) 
is decaying to two vector states ((jry), with the photon 
having no helicity zero component. Finally, to provide 
further sensitivity to B decays in the event, we admit 
two more variables. The first is the longitudinal sepa- 
ration between the decay vertex of the B meson candi- 
date and the vertex of the other charged particles in the 
event, which tends to be non-zero for signal events due 
to the long lifetime of the B meson, whereas all particles 
in continuum events usually originate from a single ver- 
tex. The second is the net flavor of the rest of the event, 
N r = N K o - Nr\, where N K o is the number 

of reconstructed K® and Nf 1 are the numbers of recon- 
structed charged particles of type i = e ± ,fi ± , or 7r^ ow 
|l2| . The discrimination power of this variable can be 
seen by noting how the final state particles (as seen by 
the detector) are produced. If the final state particle is 
produced through the decay of a B meson, then it most 
likely was generated through the weak interaction, which 
allows for flavor-changing currents. In contrast, if it is 
part of continuum u, d, and s quark-antiquark events, 
then there is no net flavor production. 

The neural network is trained on samples of simulated 
signal and continuum events. A B — > Dtt control sample 
was used to determine the systematic error in the selec- 
tion efficiency of the neural network. Here, the "bachelor 
pion" is used in place of the signal photon. This sam- 
ple is representative of signal events because the input 
variables are mainly derived from particles not associ- 
ated with the B candidate. The exception is |cosf?ff|, 
which has a distribution for pions different from that 
for photons because of the different spins of these parti- 
cles.. Instead, this variable is drawn from the signal dis- 



tribution of sin 2 0. Fig. shows the distribution of the 
neural network output for MC-simulated B° — > D~tt + 
events compared to B° — > D~ir + events reconstructed 
in on-resonance data. Also shown is the distribution 
of the neural network output for MC simulated contin- 
uum background events compared to events in the on- 
resonance sideband. These distributions show reasonable 
agreement between MC and data samples. 

The selection is optimized to achieve the best up- 
per limit because of the very small SM expectation for 
B(B° — ► 07). The quantity minimized by the optimiza- 
tion is N/e, where N is the average 90% confidence level 
upper limit on the number of reconstructed signal events 
and e is the signal efficiency. N is estimated from an en- 
semble of experiments with a given MC-simulated back- 
ground and no true signal 01 ■ We optimize the selection 
on the neural network output and on the 4> mass window 
simultaneously. The photon selection is identical to that 
optimized for the measurement of B — » K*-f 9] . The con- 
tinuum background within the signal region is estimated 
with likelihood fits to the AE* and toes distributions of 
simulated continuum events for each set of selection crite- 
ria considered during the optimization. The probability 
density functions (PDFs) used to describe the toes arid 
AE* distributions are an ARGUS threshold function 
and a first-order polynomial, respectively. The resulting 
functions are integrated over the signal region to esti- 
mate the background. The optimized selection criteria 
have a signal selection efficiency of (14.4 ± 0.1)% and a 
mean continuum background yield from MC simulation 
of 4.4 ± 0.5 events. 

The selection leads to an expectation of 0.073 ± 0.004 
BB background events within the signal region, as deter- 
mined from the MC simulation. The decays B° — > cfnr 
and B Q — > <f>r) contribute most of the events; these charm- 
less modes have not been observed, so the largest branch- 
ing fractions predicted within the SM are used to com- 
pute the expectations [15j. The remaining contributions 
come from B° — > (prj' , B° — > cf>K® and B° — ► 4>K Q L decays. 
The contributions of all other BB events, including those 
with B° — > K*°-f, were found to be negligible. 

To make a more robust estimate of the continuum 
background we use the data themselves; this eliminates 
any uncertainty related to the simulation of the qq contin- 
uum. To define functions that describe the background 
distribution we use a maximum-likelihood fit to the data 
within sideband regions I and II with the same form of the 
Toes and AE* PDFs as those used in the optimization 
procedure. The resulting functions are integrated over 
the sideband region III and the signal region to estimate 
the amount of continuum background within them. To 
validate the method the number of background expected 
in region III (mean of 79 ± 6 events) is compared to the 
number of events observed in this interval (71 events); 
these two numbers are in good agreement. To assign the 
error on the background estimate, we take the difference 
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of both PDFs evaluated with the central value of the 
fit parameters and the parameters varied by one stan- 
dard deviation and then add the differences in quadra- 
ture. The mean continuum background estimate within 
the signal region is 6.0 ± 0.9 events. 

In addition, we perform other consistency tests to en- 
sure that both the background estimate and the method 
used to acquire it are reasonable. We use a second-order 
polynomial as the PDF for AE*; this leads to a mean 
continuum background estimate of 5.6 ± 0.8 events. The 
difference between 5.6 and 6.0 (the number obtained with 
the method defined previously) is taken as a systematic 
uncertainty on the background. To check the assump- 
tion that the PDFs factorize, we compute the ARGUS 
parameter and the slope of the first-order polynomial in 
different intervals of AE* and ibes! t ne parameters are 
found to be in reasonable agreement. We perform a fit 
to Toes and AE* using regions I, II and III; this gives 
a mean expected background of 5.2 ± 0.8 events, which 
is in agreement with the expectation from using only re- 
gions I and III. The mean continuum background esti- 
mate within the signal region including the systematic 
uncertainty is 6.0 ±1.0 events. 

We observe 8 events in the signal region, which is con- 
sistent with the background estimation. The distribution 
of data events in toes and AE*, together with the pro- 
jections in these two variables, can be seen in Fig. 0] 

The fractional systematic uncertainties on e, the num- 
ber of B°B° pairs in the data set (N BO g ), and B((f> — > 
K + K~) are summarized in Table their origins are 
briefly described below. The difference in tracking ef- 
ficiency between data and simulation is evaluated with a 
sample of tracks that are well reconstructed in the SVT; 
these studies yield a relative shift of (—0.8 ± 1.3)% per 
charged track. The uncertainty related to the kaon iden- 
tification is evaluated with simulation and data samples 
of the decay D* + — > D°(— ► K~tt + )it + and its charge con- 
jugate. The photon identification efficiency is studied in 
samples of 7r° decays from t + t~ events. The uncertain- 
ties and corrections due to the ir° and 77 vetoes and shower 
isolation criteria are evaluated by embedding simulated 
photons into both data and simulation samples of BB 
events. The difference between the neural network selec- 
tion efficiency in the B° — > D~tt + data and simulation 
samples is (— 1.5±2.7)%. There is a small uncertainty re- 
lated to the statistics of the signal simulation sample used 
to calculate the efficiency. The corrected value of the ef- 
ficiency after all corrections are applied is (13.9 ± 0.7)%. 
The number of BB events used in the analysis is cal- 
culated from the change in the ratio of multihadron to 
events between data taken at the T(4S) resonance 
and at e + e~ CM energy 40MeV below the resonance; 
there is 1.1% fractional uncertainty in this measurement. 
We assume that half the number of BB events are B a B° 
events. The measured branching fraction B((f> — > K + K~) 
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FIG. 4: In (a), the t/ies - AE* distribution of data events 
after all selection criteria have been applied is shown. The 
box with the solid line indicates the signal region, while the 
dashed lines indicate the regions defined in the text. In (b) 
and (c) are the A_E* and tties projections respectively of (a). 



has a fractional uncertainty of 0.6% Q. 

Using the signal efficiency, B{(j> — > K + K ), N BO go, the 
background estimation along with the associated uncer- 
tainties, we find by the procedure of Ref. the upper 
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Systematic effect 


Correction Uncertainty (%) 


K + K~ tracking 


0.984 


2.6 


K + K~ identification 


- 


2.0 


Shower separation 


_ 


2.0 


7t°/t7 veto 


- 


1.0 


Photon detection efficiency 


0.997 


2.5 


Continuum suppression 


0.985 


2.7 


Simulation statistics 




0.1 


Overall signal efficiency 


0.966 


5.4 


BB counting 




1.1 


K+K-) 




0.6 



TABLE I: The sources of systematic uncertainty on estimat- 
ing the branching fraction. Any correction factor applied to 
the signal efficiency related to the systematic effect is also 
given. 



limit: 

B(B° -> </ry) < 8.5 x 1(T 7 , 

at the 90% C.L. In conclusion, no evidence for the decay 
B° — > 07 is observed. We set an upper limit that is 3.9 
times lower than the previously published result. 
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